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Abstract 

We discuss how to calculate Yukawa couplings in string-derived orbifold GUTs. As 
an application we investigate the quark mass hierarchies in these models. An interplay 
of different mechanisms derived from string theory leads to an interesting pattern. We 
discuss concrete examples in the context of heterotic orbifold compactifications. 
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1 Introduction 



In the standard model (SM) of particle physics the structure of Yukawa couplings of quarks 
and leptons is completely undetermined. Yukawa couplings enter the model as free pa- 
rameters. Adjusting these free parameters in an appropriate way, the theory is extremely 
successful in describing experimental results. Unfortunately it cannot explain the differences 
in the couplings and mixings. It would be appealing to explain all these input parameters 
like masses and mixings from a lower number of free parameters. 

Heterotic string theory is one of the candidates for completing the SM in the ultraviolet. 
It was shown in the last years that heterotic orbifold compactifications can reproduce the 
particle content of the SM or its minimal supersymmetric version (MSSM) [Il[2]- Yukawa 
couplings are in these models no longer free parameters but depend on the size and the 
shape of the extra dimensions and the string coupling constant. It is possible to calculate 
the Yukawa couplings from first principles and compare them, as a test of the models, with 
the experimental data. If such a model is to describe the real world, it should reproduce 
the observed quark and lepton masses, CP phases and mixings. Recently it has been shown 
that it is possible to obtain the correct value for the top Yukawa coupling in these models 
quite naturally [3]. 

Motivated by this success this study is devoted to a calculation of the Yukawa couplings for 
the other particles. It is nearly impossible to calculate the couplings in string theory exactly 
due to the complexity of the calculations. Nevertheless the ratios between the couplings can 
be computed in a simple way. We will consider mainly tree-level effects, which means that 
we will ignore quantum corrections such as those coming from the renormalization group. 

The comparison of the ratios between the different quark masses with the experimental 
data is quite exciting. We will show that the so-called benchmark models lA and IB [2] seem 
to be excluded due to their flavor structure. Maybe loop corrections induced through SUSY 
breaking |3] are able to break the obtained relations in the right way to get experimentally 
viable models. 

In section [2] we will briefly introduce the mini-landscape models before we show in section 
[3] how we calculate the ratios between different couplings. In sections H] and [5] we elaborate 
the flavor structure of the models lA and IB from [2j in more detail. Some technical details 
can be found in appendix |X] and [Bl 

2 The mini-landscape models 

In [2j several string derived models have been presented which are based on the so-called 
Zg-ii heterotic orbifold compactification. It was shown that models with the exact MSSM 
spectrum with additional heavy vector- like exotics exist. These models have appealing phe- 
nomenological features like matter parity, a large top Yukawa coupling [3] and the possibility 
to include the seesaw mechanism to create neutrino masses [5j. The flavor structure of these 
models was also investigated in [6]. We will focus here on a specific example and show that 
the structure of the models is quite restrictive in a concrete example. 

We focus on models which allow for a 6D orbifold GUT limit like the ones considered in [3] 
or [7]. Throughout this analysis we will assume that extra dimensions are compactified in 
such a way that this limit is realized. The geometry of the two extra dimensions is T^/Z2. 
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The emerging extra dimensional structure is depicted in figure [TJ The SM famihes are located 




SU(5) X U(l)^ SU(4) X SU(2)l x U(1)' 
l« 7ri?5 H 

Figure 1: A 6D orbifold GUT. The fixed points are labeled with the two quantum numbers 
(712,77-2). The first family sits at the fixed point labeled with (0, 1), the second family at the 
fixed point (0, 0). The third family and the Higgs fields live in the bulk. 

at different points in the extra dimensions (see figured]). These different location properties 
are not introduced by hand, but are a consequence of several string selection rules (see for 
example [8]). We will see that the flavor structure mainly depends on this localization. The 
flelds are further localized in the additional four extra dimensions needed to have a consistent 
string theory. The influence on the flavor structure of these extra dimensions is sub dominant 
as we will see in the next section. We can thus completely focus on the 6D GUT structure. 
A detailed discussion of the complete geometry can be found for example in [S]. 



3 Yukawa couplings in the heterotic mini-landscape 

Yukawa couplings in the mini-landscape are not only 3-point functions, but are induced by 
higher order terms due to SM singlets. Let us discuss the origin of the Yukawa couplings in 
a concrete example. Consider the up type quark couplings 

Wvukawa D Y:^^{R, s)qiu,4>i (3.1) 

where s denote standard model singlets, R represents the dependence on the extra dimensions 
and the fields are labeled according to the notation in [2]. From a bottom-up perspective 
Y^^{R, s) are free parameters which can be used to generate the observed quark masses. By 
contrast, in string theory these couplings are calculable. Matching string scattering ampli- 
tudes to the couplings in the corresponding low-energy field theory results in the schematic 
relation 

Y^{Rrs)=A^{R)-VEys (3.2) 

where A{i{R) is the string scattering amplitude. The expression "VEVs" refers to the vacuum 
expectation value of SM singlets s which induce this coupling. In general Y^^ {R, s) is a sum 
of several contributions (see for example equation (13.31) ). A{i{R) is the correlation function 
of string vertex operators on the string worldsheet which can be evaluated with the help of 
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conformal field theory (CFT). It depends on the quantum numbers of the fields and especially 
on their origin in the extra dimensions. 

Let us specifically consider the so-called benchmark model lA of [2] . The Yukawa couplings 
for the first generation of up quarks are given by 

A'u = ■ {hih^ + s?7S?8 + + + + > (3.3) 

where (p^ is the up type Higgs field and the additional fields are SM singlets (s). The 
brackets denote that we have to perform a CFT calculation to obtain a numerical value. 
The fields inside this correlator have to be understood as labels for the corresponding vertex 
operators on the string worldsheet. Because also the quarks and the Higgs field participate 
in this calculation we have included them in the above formulae. To avoid confusion we 
want to note that we use the same bracket notation for the VEV of a field as well as for 
string correlation functions. The value of the correlator f l3.3p together with the VEVs of 
the singlets is then the model prediction for the coefficient Y^^{R, s) given in equation f l3.ip . 
The singlet VEVs can be understood as expansion parameters of the superpotential. 

The calculation can be performed using CFT techniques. The final result depends only on 
the localization of the fields in the extra dimensions. The reason why only the above outlined 
couplings are allowed are several string selection rules specified in detail in [8] and [2j. 

We want to explicitly calculate the above coupling in string theory in order to explain 
or at least motivate the hierarchy between the quark masses in the SM. The techniques to 
calculate such couplings in string theory have been developed long time ago [9l[Tnj . However, 
a technical obstruction to performing such a computation results from the fact that couplings 
in the heterotic mini-landscape are generated via SM singlets and we have to deal with 
higher order couplings. There are calculational tools available also for higher order couplings 
[TTHTl] . but it is difficult to obtain numerical quantities from these general results. Excited 
twist operators which appear through oscillator modes have not been considered in the 
literature for higher order couplings so far. These couplings occur frequently in the heterotic 
mini-landscape and it seems that we have to develop tools to address these couplings as well. 

Before getting lost in a lot of technicalities it is worth to look at the structure of couplings 
first of all. We will see that it is not necessary to develop new tools, because the coupling 
structure is quite simple. As we shall see, the coupling strengths mainly depend on the 
location and on whether the fields have oscillator excitations. Table [T] shows the relevant 
data. The notation is the same as in [2]. The six extra dimensions of the Zg-n orbifold 
compactification factorize in three two dimensional tori = Tq^ x Tgy^-g^ x TgQ^-^^. Each 
torus is described by a root lattice and an orbifold point group. As already mentioned we 
consider the limit in which the torus with the SO (4) root lattice and the point group is 
large compared to the other four extra dimensions. The quantities N^. indicate if a state 
is from the oscillator sector in the corresponding torus. Here /i refers to the G2 torus, 
whereas /2 labels the SU(3) torus and /s labels the S0(4) torus. The G2 and SU(3) torus 
form four extra dimensions, whereas the SO (4) torus gives us the two additional extra 
dimensions depicted in figure [TJ No state is an oscillator in the S0(4) torus which simplifies 
the calculation as we will see later on. The quantities n2 and n'g describe the localization in 
the S0(4) torus as already discussed and describes the localization in the SU(3) torus. A 
star (*) means that the state is a bulk field in this torus. The localization in the G2 torus is 
not displayed in the table. 
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Table 1: We show a set of relevant quantum numbers for some states of the model lA. For 
the calculation of couplings the geometric origin and the oscillator numbers are highly 

important. 



We can see that not all fields are twisted and not every field is an oscillator state. As 
already remarked, computing all involved correlators in string theory exactly will be painful. 
A useful strategy to get phenomenological interesting results has been examined for example 
in [3]. Computing a quantity exactly is often quite difficult, but the relation to another 
observable is sometimes very easy to obtain. We will use a similar strategy here. We do 
not calculate the Yukawa couplings exactly, but look at the ratios between them. As we 
can see by looking at the correlators in string theory, it is irrelevant where a term comes 
from, for example from a singlet or a quark. The only important point is where the fields 
are located. The reason for this behavior is of course that in heterotic string theory the 
different particles are just excitations of a closed string in ten dimensions. Because we 
are considering the interactions between closed strings we already incorporate all essential 
features. An additional advantage of this approach is that we do not have to take care about 
the Kahler potential because it is universal for the different couplings (see also [15]). 

To be more concrete, let us take a look at another coupling. We find 

•^l^ = {(ll^2(Pl ■ sls^iS^Q [hih2 + s\jS\^ + S20S22 + ^21^2 + Ao^2-i + ^21^2'i)) ■ (3-4) 

The structure is the same as for the Y^^ coupling except for the fact that was interchanged 
by s\^. These singlets have the same quantum numbers (as can be seen in the appendix of [2]), 
but live at different fixed points (see also table [1]). This is not an accident but a consequence 
of a Z?4 symmetry [6]. The only difference between ui and U2 is, that they also live at 
different fixed points. This means, that the interchange in the singlets just compensates the 
different location of the second family. We get 

= ^My^'. (3.5) 
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Particle Twisted Oscillator N}^ N*^^ N*^^ rig n2 n'^ 
/ - -'--000 

/ - - - - 

/ / 1 - - 

/ / - 1 - 1 

/ / 1 - - 1 

Table 2: More information on the participating fields in the up quark sector. 

With (s) we denote the vacuum expectation value (VEV) of a field s. We obtain similar 
relations like that, which only depend on the VEVs of some singlet fields and not on some 
difficult string theory calculation. Using these relations we can determine the structure of 
the Yukawa couplings in dependence of the yet unknown quantity Y^^. 

4 The flavor structure of model lA 
4.1 The up quark sector 

We will consider the up quark sector in the so-called benchmark model lA in the heterotic 
mini-landscape in this section. In principle, this analysis can easily extended to all models 
in the mini-landscape. In [2] the couplings in the up quark sector were given by 

. (4.1) 

We have already looked in more detail at the couplings Y^^ and Fj^. The entry Y^^ was 
set to 1 because it already appears at the renormalizable level. In [3] it was shown that 
this coupling comes from the gauge coupling in the higher dimensional gauge theory and 
is slightly suppressed due to locahzation effects to fit in the right value. For our purposes 
we can just assume that the coupling to the top quark is of the right size and set it in the 
following to yt- 

We have collected the behavior of all necessary fields in table [T] and table [21 We start 
with a closer look at the coupling Y^^. As can be seen from the quantum numbers in table 
[T]the fields are located at different fixed points in the S0(4) torus depicted in figured] This 
different location will be encoded in the quantity /(-Re) 

Yu' = m)Ais). (4.2) 

We will look at this quantity in more detail later (see also appendix . We will see that 
the quark mass hierarchies are mainly controlled by /(-Re)- The quantity ^4(5) is no longer 
depending on the localization of the fields in the SO (4) torus, but still on the localization 
in the other extra dimensions and on the singlet VEVs. The localization in the additional 
four extra dimensions will be irrelevant for our discussion and is therefore suppressed. The 



<l2 

13 
^0 



'11 



6 



reason for that is that ^4(5) appears in every entry of and can therefore not lead to a 
difference between the quark couphngs. 

Looking at the other couphngs in detail we get 



Y 



21 



Y 



12 



Y 



22 



ACs) 



(4.3) 



where we have taken the different localization of the second family in into account. We 
have the same fields as for the Y^^ coupling despite the two quarks which are located at a 
different fixed point. For the Y^"^ coupling all fields are located at the same fixed point in 
the SO (4) torus and we get no contribution from the localization in this torus. The coupling 
is completely insensitive to the volume of this torus. 

When including also the third family things get more complicated. The third family is not 
localized and has its origin in the untwisted sector (see table [2]). In order to obtain couplings 
to localized states, extra singlet VEVs are required, such that the couplings appear at higher 
order. We obtain for example 



/lis _ jj-I oO„0„0 ^0 (U U _L_ „0 „0 I „0 „0 I „0 „0 I „o „0 I „0 „0 
~ \yi"3Y^l ''4''6''11''26 V''l"'2 "T ■='17'='18 ' ''20''22 ''21''22 ''20''23 ''21*23 



(4.4) 



Compared to Yu the singlet Sg is here playing the role of Sg whereas s^i is playing the role 
of Ml. Because M3 is from the untwisted sector it is irrelevant for our considerations. There 
is an additional term at order five in the singlets which we cannot easily relate to Y^^. We 
obtain 



Y 



13 



(s'e) (^?i)..ii 



'A) 



Y^' + B{R^rs) 



Y: 
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11 



(4.5) 



where B{Rq, s) stands for an unknown contribution appearing at order five in the singlets. 
Furthermore we have 



Y 



23 



Y 



32 



(4) {s'n)_YU 



(4.6) 



The last term occurs already at the renormalizable level and we set it to Y^^ = yt as already 
discussed. Putting everything together we get the matrix 



Y, 



f{R6)Ai~s) 



4) 



f{R6)Ai~s) 



'-^f{Re)Ai~s) + B{R,,S)\ 



(4) 

(.0) f{Re)ACs) 



ACs) 



4) {4,) 
(4) 



f{R6)ACs) 



yt 



(4.7) 



This matrix contains several parameters. The singlet VEVs are, in principle, calculable as 
they are determined by F- and D-fiatness constraints. However, in practice it is quite difficult 
to compute them. In what follows, we will see that one can nevertheless gain non-trivial 
information since we can understand the behavior of /(-Re) from a string theory calculation. 
Another important information is that the diagonal entries of Yu are completely independent 
of the VEVs. Let us now discuss the quantities f{Re) and B{Rq, s) in somewhat more detail. 
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Quark Mass 
up quark (ui) 1.5 — 3.3 MeV 

charm quark (^2) 1-27 ( ] GeV 



^-0.11^ 

top quark (wg) 171.3 ± 1.1 ±1.2 GeV 
Table 3: Up quark masses from [T6] . 

The expression /(-Re) reflects the different locahzation of the flrst family in the SO (4) 
torus. It thus depends only on the moduli dependent part of the amplitude in this torus. 
We have four twisted flelds and no oscillator. Looking at the location of the flelds we can 
see that for Y^"^ all flelds live at the same flxed point in this torus. We conclude that there is 
no suppression due to instanton effects (see for example [9j). For Fj^ the situation changes. 
We get differently localized flelds and thus an instanton suppression. As a flrst result we can 
assume that Y^'^ > Y^^ and therefore /(-Re) ^ 1- 

Whenever two flelds from different flxed points interact the coupling is suppressed by the 
distance between the flelds. If only flelds from the same orbifold flxed point or from the bulk 
interact the coupling is unsuppressed. 

A detailed look at a four twist calculation (as for example performed in [9]) tells us that 
the only difference between Y^^ and Y^"^ is in the classical part of the amplitude. This is 
the instanton part which goes roughly like e'^^s where vr-Re is the distance between the two 
flxed points labeled with n'2 in the SO (4) torus (of course this distance is encoded in the T 
modulus). We get the following suppression (see appendix |A] for a detailed calculation) 

/(i?6)=e-^. (4.8) 

The quantity B{Rq,s) is quite difficult to discuss. It is a coupling involving six oscillators. 
It is also the only coupling depending on 5^3. B{Rq,s) appears only in one entry of the 
Yukawa matrix and thus it affects the result only marginaly. 

Using this ingredients we can try to obtain quantitative results. The running of the 
couphngs is neglected as it is beyond the accuracy of our approach (see also [5]). As a 
flrst test, the eigenvalues of Yu should give us the right ratio of the different up type quark 
masses. We start with flxed VEVs, all of the same size (s°) = 0.1. As usual we can change 
to mass eigenstates to relate the couplings in Yu to the quark masses. The diagonal entries 
of the rotated matrix have to be proportional to the quark masses. The experimental values 
are given in table O We get ratios between the quark masses like 500 between the flrst two 
generations and 100 between the second and the third generation 

m,, ! 1 rrif, \ 1 , ^ 

— ^ , — ^ . 4.9 

rric 500 100 ^ ^ 

The impact of the VEVs and B{Rq,s) is negligible and we ignore it. We therefore set 
B{Rq, s) to zero. 

The dependence on the remaining quantities is approximately given by 

^^/(i?6), ^^A(5). (4.10) 

rric irtt 
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Quark Mass 
down quark {di) 3.5 — 6 MeV 

strange quark (^2) 105 (~'~^^] MeV 



+0.17 

bottom quark (^3 — ^4) 4.2 ( ' ) GeV 



35 

).r 

0.07 

Table 4: Down quark masses from [TB]. The labels for the quark masses are the same as 

ind. 



If we choose /(-Re) = 5^ 'we can get the right behavior with A{s) = The off-diagonal 
terms are now quite small and we obtain the eigenvalues m ~ (0.00002,0.01, 1). Given the 
experimental values (displayed in table [3]) these values are of the right size. 

The small value of /(-Re) can be explained by the instanton suppression in the S0(4) 
torus. Using our result from equation fl4.8l) we obtain Rq ^ 2.8. As also remarked in [3] 
we expect to relate the distance -R5 between the two fixed points labeled by n2 to the GUT 
scale, namely Mqut = As a starting point we use -R5 = 50, which gives us Mqut = 
in Planck units. With these values we find an anisotropy of a factor of = 18 to get the 
right value for -Re. Given |3] such a value seems to be excluded. Nevertheless it is very 
encouraging that we get a value in the correct region. As outlined in section 14.21 the model 
suffers from more serious problems, therefore we do not address this point here in more 
detail. Another uncertainty is the value for A{s) which may or may not be realistic. 

We conclude that even if we have a lot of unknowns, we can use the strong correlations 
between the couplings due to the different selection rules to make predictions. In contrast 
to the general analysis of [6] we find that the models are quite predictive. It is appealing 
that we have to consider in principle only a 6D GUT with differently localized families to 
explain the hierarchy for the up quark masses. We have obtained a very intuitive geometrical 
explanation and no miraculous string features at work. Taking the result of [3] into account 
the scenario is very constraint. We have one free parameter, namely -Re, to explain the 
ratio between the gauge coupling and the top Yukawa coupling and simultaneously the ratio 
between the up and the charm quark mass. 

We proceed by looking at the remaining quark masses in the down quark sector. 



4.2 The down quark sector 

After the success of relating the quark masses in the up quark sector to the geometry of the 
extra dimensions, we will also try to explain the hierarchy for the down quark masses on 
the same footing. As we will see one cannot fit the up quark and the down quark masses 
simultaneously. The experimental data are displayed in table IH We find a mass ratio of 
19 between the down and the strange quark mass [17J and a ratio around 40 between the 
strange and the bottom quark mass 

rrid I I nis \ I 

— ^ — , — ^ — . (4.11) 

rus 19 mb 40 ^ ^ 

We take into account couplings up to order nine in standard model singlets which correspond 
to 12-point couplings. Because we consider also couplings at higher order in singlets our 
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starting point is slightly different from [2] and we start with 



§5 



Yd= \ s-^ (4.12) 



which is a more accurate starting point. Repeating the steps of the discussion in the up 
quark sector we end up with a matrix like 

fg,{Re)C{s) fmD{~s) \ 
Yd = f{Re)D{~s) g2{R^)C{s) . (4.13) 
V X(i?6,5) U{s)) 

We have a lot of parameters, from a naive perspective it should be easily possible to fit the 
right values for all quark masses. 

For the coupling between and the massless combination ^3 — we find a coupling at 
order eight which is unsuppressed in the SO (4) torus because we find a coupling where all 
fields are localized at the same fixed point. Looking at the eigenvalues of the matrix we can 
see that the coupling U{s) exactly gives us one eigenvalue. As outlined above the coupling is 
not related to R^ or R^ and we can treat it as a free parameter to fit the bottom quark mass. 
The situation is similarly as for the top quark coupling, the third generation decouples from 
the first two generations. 

For the other two quark masses the situation is more complicated. We find that X{Rq^ s) 
plays no role for the eigenvalues of Y^ and is thus irrelevant for our discussion. The off- 
diagonal couplings f{RQ)D{s) arise already at order five in standard model singlets and are 
suppressed in the SO (4) torus. We find the same factor as in the up quark sector, namely 
/(-Re) for the suppression. /(-Re) is the same as in the up quark sector, we have again four 
twisted fields in the SO (4) torus. Again two fields are localized at the fixed point labeled 
by n2 = 0, whereas two of them are localized at the fixed point n2 = 1. To be concrete, we 
obtain 

= {q,d2ct>i ■ slslsl, {hh, + s',,s% + s%sl, + sl.sl, + s^.sl, + s^.s',,) > + Oi~s') (4.14) 
and thus 

Inserting already the Higgs VEVs we can relate quantities from the up quark sector to the 
down quark sector. We have seen that Y^^ = /(-Re)A(s) together with tan/3 gives the up 
quark mass. We thus can see that Yj^ is of the right order to give us the down quark 
mass because the singlet VEVs are supposed to be of the same order of magnitude (see for 
example [IB])- Of course, in detail the situation is more constrained, because tan/3 = 1 is 
excluded by experiments. An allowed value of about tan (3 > 2 [19] yields a too small down 
quark mass if the VEVs are of the same size. Nevertheless, this can be overcome due to an 
according VEV assignment. If this can be realistic or not is an open question. 

We further note that we get contributions from higher order terms. These corrections do 
not affect the statement Yj"^ = Yj^ which is correct also up to order nine in the standard 
model singlets. The reason for this behavior is the -D4 symmetry |6j. 
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The diagonal couplings gi{RQ)C{s) and g2{R6)C{s) arise at order nine in standard model 
singlets. Their calculation is challenging. We have (for the quantum numbers of the addi- 
tional fields see [2]) 

= 9,{Re)C{~s), Yr = ^^Yj' (4.16) 

with 

= ((lldl(pl ■ S4 (s?3)^ 526^30 (XlX2/i2/i4 + XzXih2K + XiX4h2he + X2X3h2h6 

(4.17) 

+h2hs{h^f + h2h'i{hQf + h2hJi^hQ) 

We find a similar behavior as in the up quark sector. The couplings for the first and second 
generation are not the same, but differ due to the different localization of the second family. 
We introduce the factors gi{R&) and (72 (-Re) which encode the localization behavior. We find 
for these couplings eight fields twisted in the SO (4) torus. There are always fields localized 
at different fixed points. We conclude that despite the fact that the coupling is already 
suppressed by a high power of singlet VEVs, the coupling is further suppressed due to the 
localization in the S0(4) torus (see also appendix [B]). Compared to the off-diagonal terms, 
these couplings are therefore rather small. We have to introduce two different factors because 
both couplings are suppressed. In the up quark case only the Y^^ coupling was suppressed. 

We want to look at the eigenvalues to compare them to the experimental quark masses like 
in the up quark case. X{Rq, s) is irrelevant for the eigenvalues due to the zeros in the matrix 
Yd- Because the off-diagonal entries are identical and compared to the diagonal couplings 
quite large we get two almost degenerate eigenvalues. This is in contrast to data. To be more 
concrete we can estimate the difference between the couplings. If we assume rather large 
VEVs (5) = 0.3 we get a suppression of (s)"^ ~ between the diagonal and the off-diagonal 
couplings because the diagonal couplings arise at order nine and the off-diagonal couplings 
at order five in standard model singlets 

C(5) ^ ^D{s). (4.18) 

We should point out that for smaller VEVs this ratio is even larger. If we use the results of 
appendix [B] we obtain 

fJ^J{R,)D{s) f{R,)D{s) \ 
Yd ^ f{Re)D{~s) ^J{R,)D{s) . (4.19) 
V X(i?6,5) U{s)j 

We can tune D{s) to obtain the right ratio between the bottom and the strange quark mass 
regardless of Rq, but it is impossible to get a big ratio between the down and the strange 
quark mass. The quark mass ratios are 

m,^62 rn^^VmfJMm 

nis ^ 63' mfe ^ 125 U{s) ' ^ ' ' 

The mass eigenvalues for rud and are almost degenerate for all possible values of the free 
parameters. 
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We conclude that it is impossible to obtain a quark mass hierarchy as expected from 
experimental data with our setup. 

It was first proposed in [4j to generate Yukawa couplings via loop effects through SUSY 
breaking soft terms. This can be a mechanism to cure the problem in the down quark 
sector. This mechanism leads also naturally to flavor changing neutral currents (FCNCs) 
which have to be compatible with experimental bounds (see for example [20] in the string 
theory context). If it is possible to fulfill these bounds in this class of models remains an 
open task for future work. Maybe there are also other mechanism which can give rise to 
essential corrections to the down quark masses in the right way. 



5 A comment on model IB 

In [2] another vacuum configuration called model IB was considered. This model has an 
unbroken D4 family symmetry for the first two generations. The difference to model lA 
is that the symmetry is not broken by the VEV assignment. The consequence is that all 
couplings for the first two generations are equal. As we have seen in the last section for the 
model lA, such a behavior is in conflict with observations. It seems to be impossible for 
such configurations to get a big ratio between the masses as required by experiment without 
radiative quark mass generation already mentioned in the last section. We conclude that 
these configurations are rather unrealistic by construction. 



6 Summary 

We have developed tools to study the quark mass pattern in heterotic orbifold compactifica- 
tions. Specifically we have studied the flavor structure in two examples of the mini-landscape, 
namely the so-called benchmark models lA and IB. We have seen that without any difficult 
string theory calculation we can obtain several interesting results because the couplings enjoy 
relations. We have linked the quark masses to the extra dimensions like already in [3] for 
the top quark coupling. We have also found some relations between the up quark and the 
down quark sector depending on tan/3. Our study shows that it is possible to obtain simple 
relations for the Yukawa couplings even in complex string-derived setups. 

For the concrete models studied in our analysis we found that it is impossible to get 
the right masses in the down quark sector without including large radiative contributions 
through SUSY breaking. We do not know if this feature is universal for all models in the 
mini-landscape. What we can state is that it seems to be quite easy to disfavor models 
from the mini-landscape. Because we have not succeeded in obtaining the right quark mass 
pattern, we have not looked at the mixing, namely the CKM matrix or the mass pattern in 
the lepton sector. With the techniques outlined in this paper this is in principle also possible. 

In [21] and [22] powerful tools have been developed to study the discrete symmetries 
between the string selection rules. Maybe it is possible to find better string theory models 
with these techniques. 
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A AT-point couplings for four %2 twists 

To obtain the quantity /{Rq) we need to calculate a A^-point coupling with four twists. 
The calculation of a 4-point coupling with four twists has been calculated long time 
ago [HlllO]. Let us review here the result. There are two different contributions, one from 
the quantum part and one from the classical part. The first is insensitive to the localization 
of the fields, whereas the latter one is sensitive. Thus the quantity /(i^e) is only influenced 
by the classical part. The result for the classical part is [H] 



where we have 

2Fi(i,i,l,x) K{x) 

and ej labels the fixed point of the twist fields. With K{x) we mean the complete elliptic 
integral of the first kind, x can be interpreted as the position of one of the vertex operators 
on the sphere. It cannot be fixed by PSL(2, C) invariance and has to be weighted by an 
integral in the end. That means we have one free complex modulus of the Riemann surface 
describing the string interaction. We further have = for the fixed point (?t,2, ^1-2) = (0, 0) 
and €4 = 1 for the fixed point (n2, n'^) = (0, 1). 

If all fields live at the same fixed point eo = ei = we get an unsuppressed amplitude 
because the exponential breaks down. If the fields are localized at different fixed points, for 
example eo = and ei = 1 we get a suppression going with the distance squared R^. The 
terms for 7^ are further suppressed and come from states winding also around the torus. 
We can neglect their contribution for simplicity and obtain 

all fields live at the same fixed point 
Zci{x) ^ { .nl . (A.3) 

e fields live at different fixed points 

This is the exact classical contribution in the SO (4) torus, because all other fields are un- 
twisted in the SO (4) torus. We do not know the contribution in the other extra dimensions, 
namely the SU(3) and the G2 torus in detail. We know that the additional contribution is 
depending on x. Because we have to integrate over x in the end, we have to approximate 
the result. 

Let us focus now on the A^-point coupling in more detail. We have a A^-point coupling 
and therefore A^ vertex operator positions on the worldsheet. We can fix three of them by 
PSL(2, C) invariance. Over the remaining A^ — 3 positions we should integrate. At least one 
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Re(x) 

-1.0 

Figure 2: The function Sc\{x) for the values eo = 0, ei = 1 and Rq = 2. 



of these vertex operators is twisted in the SO (4) torus and we called his position x (in fact it 
is also possible that all vertex operators in the S0(4) torus are undetermined by PSL(2, C) 
invariance, but we will not consider this case here). The other vertex operators are untwisted 
in the SO (4) torus and we do not know their contribution exactly, but they depend on all 
— 3 variables. Their contribution is the same, regardless of their localization in the SO (4) 
torus. The complete contribution coming from this localization enters via equation f lA.3|) . If 
we call the unknown part from the SU(3) and G2 torus and the contributions not affected 
by the localization in the S0(4) torus U{x, X2, ■ ■ ■ xn^) we get 

J (fx J (fx2--- J (fxN-3Zcl{x)U{x, X2, ■ ■ ■ Xn-s) (A. 4) 

for the complete amplitude. We can approximate the function Zci{x) by setting x = \ which 
is motivated by the shape of the function (see figure [2]). We obtain r(x) = i for this choice. 
We get 

(fx I <fx2--- / d?XN-'iU{x,X2., ■ ■ ■ xn-z) 



if all fields live at the same fixed point and the coupling is unsuppressed. In the case where 
the fields are localized at different fixed points we get 



e 4 Jdxjdx2---Jd xn-3U{x, X2, . . . xns)- (A. 6) 

We recognize that the integrals are now the same. We conclude that the required ratio 
between an suppressed and an unsuppressed coupling is 

j,^^ ^ _ e 4^ / d'^x J d'^X2 ...J d^XN-sU (x, X2, . . . x^s) _ 

/ d'^X J d'^X2 ■ . . J d'^XN-3U{x, X2, . ■ ■ Xn-s) 

in a first approximation. The uncertainty of this approach is of course substantial, because 
the function U{x,X2, ■ ■ ■ XN-3) is unknown. 
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B AT-point couplings for eight %2 twists 



The functions gi{Re) and g2{R&) can be determined by looking at A^-point couplings with 
eight ^2 twists. The computation of this amplitude is more complicated than the one 
reviewed in appendix |Al General 'L2 twist amplitudes have been considered in several papers 
[T2 | 1231427] ■ also because they are of interest in the context of spin operators. The beautiful 
mathematics of hyperelliptic Riemann surfaces needed for these computations can be found 
for example in [2S]. The factors Qi^Rq) and (-Re) are not only given by the classical part of 
the amplitude, which is now depending on at least five free moduli of the Riemann surface 
Xi, . . .x^. The reason for this is that different fields appear for different couplings and the 
quantum part is depending on all of these fields. Because our approach contains already a 
big uncertainty we assume 

^qu = Oil) (B.l) 

and neglect the quantum part for the following discussion. This can be justified due to the 
discussion in [29] in the context of magnetized brane models and will be discussed in more 
detail for heterotic orbifolds elsewhere. 
The classical part can be written as 

Z,i = J]e-^='(^i'-'^5l (B.2) 

The dependence on Rq is the same as for f{R%). We find 

gi{R,) = e-"»^i (B.3) 

and thus again an exponential suppression of the order of f{R%)- The constants > are 
depending on the period matrix of the Riemann surface and on the integral approximation. 
The couphng (71 (-Re) is dominated by couplings where six twist fields are localized at one fixed 
point and two twist fields at the other fixed point. We thus conclude that gi{Re) ~ /(-Re) 
(see for example [I2]). For the coupling (72 (-Re) the dominant contribution originates from 
the same coupling as for the giiR^) term. There is an additional contribution where four 
twist fields sit at one fixed point and four twist fields at the other fixed point. Nevertheless, 
this contribution is stronger suppressed [12]. We thus conclude that (72 (-Re) ~ fl'i(-Re)- There 
is only a mild difference between (72 (-Re) and gi{Re) which we will neglect. For the numerical 
discussion of couplings we assume 

g^{Re)^g2m- fiRe). (B.4) 
The error of this approximation is larger as for the calculation of /(-Re) discussed in appendix 
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